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We present a model for an inhomogeneous and anisotropic early universe filled with a 
nonlinear electromagnetic field of Born-Infeld (Bf) type. The effects of the BI field are 
compared with the linear case (Maxwell). Since the curvature invaria nts are well behaved 
then we conjecture that our model does not present an initial big bang singularity. The 
existence of the BI field modifies the curvature invariants at t = as well as sets bounds 
on the amplitude of the conformal metric function. 

1. Introduction 

The advantage in exploring a cosmological model with nonlinear electrodynamics 
resides in approaching with a classical model but taking into account some quantical 
features. It is well known that nonlinear electrodynamics (NLE) is the best classical 
solution to the self-energy problem of charged particles. The aim of nonlinear elec- 
trodynamics of Born and Infeld (BI) § was to establish a model of a finite classical 
field theory without divergences. However, NLE can also be a model of quantum 
electrodynamics in the classical limit of high occupation numbers. As demonstrated 
by Stehle and De Baryshe H, theories with Lagrangians similar to the Heisenberg- 
Euler effective Lagrangians (the weak field expansion of the Born-Infeld theory is 
of this form) are more accurate classical approximations of QED than Maxwell's 
theory in the case of fields with high intensities at a fixed frequency. This fact 
makes NLE particularly interesting in general relativity, specialy in cosmological 
theories, as a simple classical model to explain vacuum polarization processes- a 
possible influence on the mechanism of the evolution of the early universe. The 
processes called scattering of light by light can be expressed in terms of a kind of 
electric and magnetic permeability of vacuum. From this point of view, the exact 
solutions of the Einstein-Born-Infeld (EBI) equations are worth to study since they 
may indicate the physical relevance of nonlinear effects in strong gravitational and 
strong electromagnetic fields. 

The Born-Infeld theory is self-consistent and satisfies all natural requirements. 
The BI Lagrangian depends on a constant 6, that is a maximum field strenght, and 
also depends on the invariants of the field. Recently BI has recover interest in the 
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context of string theory, since the low energy action Lagrangian corresponding to 
the objects called brane^can be written as a term that is BI plus a Weiss-Zumino 
or Chern-Simmons term u. Furthermore, the BI strings are solitons which represent 
strings ending on branes. 

In relation to exact solutions of the EBI equations we can mention that there are 
exact solutions for spaces possessing one timelike and pne spacelike Killing vector 
fields. They correspoiid to spherical static spacetimes □ and also to axisymmetric 
stationary spacetimes u. 

In the literature cosmological models have been proposed for the source-free 
Einstein-Maxwell equations i, 0. Charach and Malin i presented a cosmological 
model of Einstein-Rosen with gravitational, electromagnetic and scalar waves, this 
work generalized the electromagnetic Gowdy model. In all of these models the mat- 
ter fields cause an evolution significantly different from that of the vacuum model. 
The interest is to get insight on the coupled electromagnetic and gravitational radi- 
ation at early epochs of the universe. Since at that epoch the fields must had been 
very strong, the nonlinear electromagnetic effects could be of importance. These 
fields could originate from processes near the Planck era (t ^ 10~'^^s). 

The spacetime we address in this work possesses two Killing vectors, both space- 
like. It is not immediate to conclude that this spacetime is also compatible with a 
BI source. For instance it does not admit a perfect fluid source together with BI, 
which is the case in the stationary axisymmetric spacetime. This is the first atempt, 
as far as we know, of exploring BI early universes. Our aim is to determine to what 
extent the BI field can prevent or modify the initial big bang singularity. Spatially 
inhomogeneous cosmological models with perfect fiuid which do not xiriginate in an 
initial singularity have been previously found by SenovillaQ (see also U). For the 
solutions derived in the present work the initial big bang singularity appears only as 
a coordinate singularity, since the invariants do not diverge at the origin. It turns 
out that the existence of the BI field set bounds to the amplitude of the conformal 
metric function and also affects quantitatively the curvature invariants at t = 0. In 
this work we do not treat the possible transition of our model to a homogeneous, 
isotropic one. In Sec. 2 the Einstein-Born-Infeld equations are established. In Sec. 
3 the solution is presented and analyzed in the linear limit as well as the asymp- 
totic behavior. In Sec. 3.4 we remark the coordinate transformations relating the 
spacetime studied and the Einstein-Rosen metrics. Concluding remarks are given 
at the end. 

2. Einstein-Born-Infeld Equations 

The line element considered corresponds to a G2 spacetime, given by 

ds^ = _ ^ + G[e-^dy2 + e^(dx + mdy f]}, (1) 

ip'^ h s 

where x, y are ignorable coordinates (i.e. the spacetime has two spacelike Killing 
vectors, dx,dy), the metric functions depend on z and t: ip — ip{t), h = h(z), 
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s = s{t). G{z,t) and e^{z,t) are separable functions of z and t. The metrie (|l|) 
eorresponds to a space with gravitational waves propagating in the z direction, 
with their wavefronts on the plane surfaces spanned by the Killing vectors. The 
function W describes what is called the "+" polarization of gravitational waves. 
G(z, t) describes the transverse scale expansion created by the energy density of 
the waves. The function to(z) = niQ + 2lz is playing the role of unpolarizer for the 
propagating gravitational waves that constitute the background. Making I — > the 
metric recover the form of a space of linearly polarized gravitational waves (diagonal 
metric) . 

To interpret line element (|l]) as cosmological, we refer to the Einstein-Rosen 
metric Il3, Eq. (^), Sec. (3.3). The local behavior of these spacetimes is defined 
by the gradient of the metric function G, i.e. w hich can be spacelike, null 
or timelike. The globally spacelike case corresponds to cylindrical spacetimes. The 
case of being globally null corresponds to the plane symmetric waves. The other 
case of G,p being globally timelike as well as cases in which the sign of the gradient 
can vary from point to point are used to describe cosmological models and also 
colliding gravitational waves 0. It must also be remarked that due to the presence 
of an Abelian subgroup G2 in the Bianchi models of tvpes TVII, the Einstein- 
Rosen metrics include these models as particular cases t3. Also the rotationally 
invariant Bianchi models of types VIII and IX belong to the generalized Einstein- 
Rosen spacetimes 0. 

We note that in the case of Einstein-Rosen spacetimes, the interpretation as 
gravitational waves propagating in the z direction arises directly from the equation 
satisfied by G, 

G,tt-G.^,, = 8TTGef{T^-T), (2) 

For vacuum or when = T, G satisfies the usual wave equation. In our case, 
the B-I field is not traceless and then G is some propagating wave with a source 
(the B-I field). 

Turning back to the spacetime (0), it is filled with a BI field with energy- 
momentum tensor given by: 

= 27^.pP^„P", - g^42Pn,p + Qn,Q - H], (3) 

where H is the Born-Infeld structural function and "hCx '■— dH/dx. P and Q are 
the invariants associated with the Born-Infeld field, given by 

P^^Pf^.P'"', Q^^P^^P'^'', (4) 

P^"" denotes the dual tensor of P^'',defined by P^"' := ie^^^^Pc/j- The antisym- 
metric tensor P^^ is the generalization of the electromagnetic tensor Ff^^. The 
structural function is constrained to satisfy some physical requirements: (i) the 
correspondence to the linear theory [H(P, Q) = P + 0{P^, Q^)]] (ii) the parity con- 
servation [H(P, Q) — HiP, —Q)]; (ui) the positive definiteness of the energy density 
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{TL,p > 0) and the requirement of the timehke nature of the energy flux vector 
{PTi p + QTi^Q — 7i > 0). Such structural function for the BI field is given by 



n = b^- V&^-2fe2p + Q^ (5) 

where 6 is a maximum field strcnght of the BI field. In the linear limit, which is 
obtained by taking & — > cxd, then Ji — P and P^i, ~ i^^jy, coinciding with the Maxwell 
electromagnetism. P^^ and -F^j/ are related through the material or constitutive 
equations: 

Ff,,=n,pP^.+KQP^,, (6) 
We have worked out the problem in the null tetrad formalism 0. The function 
H in Eq. d) involves a square root that complicates the solving of the system of 
the Einstein-Born-Infeld equations. To surmount this difficulty, first of all we have 
aligned the two non-null components of P^^ along the two principal directions of the 
tetrad vectors and (the metric is type D); second, we parametrize the BI field 
in terms of the electric and magnetic fields: P34 = D and P12 = iH where D and H 
stand for the usual electric displacement and magnetic field intensity, respectively. 
Also, we introduce the function 1/ = viD, H), 

e-P[2H = ^T3^, (7) 

If the previous Eq. (Q) has to make sense, it demands that b > H. From the 
structural Eqs. we obtain the expressions that clearly manifest the vacuum 
polarization, in this case depending on the function z/(z, t): 

D = e^E, B = e'^H, (8) 

Besides the Einstein field equations we must consider the BI electrodynamical 
equations that amount to the closure condition of the electromagnetic two form, uj, 

doj = d{{Fi2 + P3M A + (P34 + Pi2)e3 A e^} = 0, (9) 

This formalism of the Born-Infeld field was introduced by Plebafiski et a]0. We 
have assumed that h{z) = a + (3z + ez^, where e can take the values —1,0, 1. For 
each value of e we obtain different curvatures in the z direction. For the metric (|l|) 
and with the alignment chosen the electromagnetic equations are 

(E + iH),, = 0, (10) 

{\n{D + iB))^t = -2ile-'' + 2{lnip)^t, (H) 

Eq.(9) tells us that neither E nor H depend on z. Eq.(lO) and Eq. (^5|), quoted 
below, lead to infer that neither D nor B depend on z. Therefore, for this spacetime, 
the BI is a spatially homogeneous field. The Einstein equations for the BI energy- 
momentum tensor (^ amount to 

(^,„ + ZV = 0, (12) 

2 

^s,„ - 2(/.^,tS,t + s(3(^2 + 3ZV') = 2b^{e-''-l), (13) 
ip(p^tS.t + s{-3^^t + P^^) + eif^ = -2b^(e''-l) (14) 
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It is a system of equations for the functions f(t), v{t) and s(i), while the rest 
of the metric functions are = {s/h)^ and G{z,t) = {hs)^. Eq. (11) for (p{t) 
can be solved immediately, giving (p{t) — A cos{lt) + B sin(^t), with A, B constants. 
Eqs. (12)-(13) can be decoupled and it is obtained the solution for v{t) in terms of 
ip{t) as well as Eq. ( p^ fo r s{t): 

v{t)^\Hl-^{tf), (15) 



ifip^ts.t + s(-3¥'' + ; V) + - -262(V(l-¥>4) - 1), (16) 



In order that v{t), given in Eq. (15), be a real function, Lp{t) — Acos,{lt) 



Bsm{lt) must be such that \ip\'^ < 1. This means that the existence of the BI field 
restrains the conformal metric function ip(t) in its amplitude because it bounds the 
values of the constants A and B. Demanding that 

A<1 if B = 0; 
B <1 if A = 0; 
A^ + B^ <l if A^O, B^O, (17) 

fulfils the requirement < 1. Putting the NLE functions, D{t) and B{t), Eqs. 
(H), in terms of (p{t) 

D{t) = ip^cos[2l j {I- ip'^Y^dt], 

B{t) = -ip^sm[2l J {l-(p^rUt], (18) 
with the corresponding coordinate components of the electromagnetic field F^i^, 

F,y{t) = -e''sm[2l J {l-^^)--2dtl (19) 

The previous expressions show that the BI field modulates the amplitude and fre- 
quency of the electromagnetic field. The fields B(t) and D{t) are shown in Fig. 1. 
They are oscillating functions depending on the constant b through the function e'^. 
Since the metric is of Petrov type D, the only nonvanishing Weyl scalar is given 

by 

i^2 = ^'['-f~^~2Ps + ^^s,], (20) 

This scalar is not real and therefore the magnetic part of the Weyl tensor does 
not vanish. We also note that it is homogeneous, i. e. it depends only on time. 
The invariants are given by / = 3ip2^ J — ~4'2j from their expressions it is clear 
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that they share (if they exist) the singularities of "02 • On its turn ip2 depends on 
the behaviour of the functions ip{t) and s{t). ip{t) is a periodic function while s(t), 
governed by Eq. ( p^ ) , carries much of NLEBI information. In the next section we 
analyze the behavior of s{t). 

3. The BI solution and the Hnear case 

The solution to Eq. ( [l^ ) for s{t), in terms of (p is given by 

.w^.v.w-/ ^^''^:""-''' .^.), (21) 

where ci is an integration constant related to initial or boundary conditions for 
vacuum. It is interesting now to analyze graphics of the function s{t). Fig. 2 
shows s{t) for several values of the BI constant b. For values of 6 > 1 the function 
s{t) becomes positive for all the range and reduces its oscillations; for larger & t he 
maximums become greater. The transverse scale expansion G{z, t) ~ ^Jh(z)s{t) is 
an inhomogeneous function expanding in the z direction, it is shown in Fig. 3. 

Fig. 4 shows the plots of the real and imaginary parts of '02; they give us 
qualitative information with respect to the singularities of the invariants / and J . 
Both Re'02 and Im02 are continuous functions that do not present infinities at all. 

In Fig. 5, the function s{t) displays different behaviours at i = depending if 
(f = Bsmilt) or (fi = Acos{lt). In the former case s(0) = and consequently there 
appears a coordinate singularity at t = 0; however the corresponding Weyl scalar 
is not singular there. We can see from Fig. 5 that this class of solutions admits 
both behaviors at t = 0, with coordinate singularity or without it. To a linear 
combination cp = Acos{lt) + Bs'm{lt) corresponds s(0) 7^ 0. 

3.1. Linear limit 

In order to compare the effects of the BI field with the Maxwell case, we compute 
the linear limit. The Maxwell electrodynamics is recovered when b —f 00. The BI 
Eqs. (9) and ([ll| ) in the linear limit correspond to 

[ln(7^ + iB)],t-2(ln(^),t =0, (22) 

whose solution is 

D = E = C cos{2 Alt) , B = H = C sm{2 Alt) , (23) 

where A and C are constants. In terms of the coordinate components of the elec- 
tromagnetic field the expressions are 

Ft^ = Ccos{2lt), F,y ^ ~Csm{2lt), (24) 

This electromagnetic field is homogeneous of intensity — + . In accor- 
dance with the parametrization (|^), for b 00, Eqs. (12)-(13) become 
2 
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^ip^ts,t + s{~3ip'^t + l'^^^) + ^^^ = -C^ (25) 
Solving for s{t) in terms of (p(t) we have 

sit) = ip^ipAci - [ ^ t7 ^^]' (26) 



expression that also can be obtained from Eq.(21) taking the limit 6 — s- cx), 



lim 26^(^1 -ip-i- 1) ^ lim 262(6"^ - 1) = C^- 

b — *cxj b — >oo 

Fig. 6 shows the comparison between s(t) corresponding to the linear electromag- 
netic field, Eq. @), and s{t) of the BI field, Eq. (|2l|). The presence of BI field 
smooths the function s{t) respect of the linear version. Making C = we get the 
vacuum limit of the spacetime (|l|), by this meaning that only gravitational waves 
remain. 

Comparing Eq.(|T]) for s{t) corresponding to the BI field with Eq.(p6|) for the 
linear case we see that the first term in the integral is directly associated with the 
NLEBI effects. To show the dependence on b we write s{t) as 

sit)=^^^,[c,- I ^^^^^f dt], (27) 
J V,tV 



where f{h) = 2&^(e'^ ~ 1) = 26^{y p^'-^l — 1}, in accordance with the requirement 
in Eq. (|7|),the range of 6 is 6 > iJ. If one plots f{b) vs. b, it can be seen that the 
NLEBI effects are sensible for values of b near H , while for larger &'s, generically it 
tends to the constant linear Maxwell field = i?^ + H"^ . 

3.2. Behavior at early times 

For t << 1 it is reasonable to approximate </?(i) as 

ip{t)= A cos{lt) + B sm{lt) K A + Bit, (28) 

Since f = const, implies s{t) — 0, we take tp k. t along with the restrictions (|ll 
imposed by the BI field. The constant I is suposed not to be large since we are in 
a low frequency regime. To solve Eq. (^l|) we also approximate \/l — Lp'^ « 1 — 
obtaining 

s{t)Ket^ + cit^ + b^[t^ + 0{t^)], (29) 

predominance of one term over the others depend on the relative values between 
ci and b. The term in brackets represents the BI contribution at early times. It is 
expected b to be large if the electromagnetic field is strong, then the BI contributi 
on can be of importance. The first term corresponds to the contribution of the space 
curvature, determined by the value of e, (1, —1, 0). The Weyl scalar V'2 approaches 
t = as 
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The contribution of the BI field goes as and higher orders in t. In the magnetic 
part of '02 the contribution is not important at very early times since comes from 
the term with e = —1, 0, 1 and goes as t^. 

We now mention the case I = 0, which corresponds to a spacetime of propagating 
linearly polarized gravitational plane waves (diagonal metric). Eq. (11) for ip(t) 
become (p^tt = 0, with solution Aot + Bq, Ao,Bo being constants. The refore, the 
behavior of (p in the case I — resembles the one for t << 1. Consequently, we 
assert that this spacetime approaches the origin in time as propagating linearly 
polarized waves. 

3.3. Behavior at t >> 1 

To complete the previous analysis, the question arises as how does the model 
evolve at < >> 1, or which is the asymptotic behavior of the solution at t >> 1, i. 
e. for times when the universe is already a causally connected one. 

The dynamics of our model is driven by the function (p{t). For t >> 1 both ip 
and ip^t become rapidly oscillating periodic functions that can be approximated by 
a constant. With this, the metric function s{t) (Eq. ( |2l| ) becomes a linear function 
on t. Absorbing constants and transforming ^/Ki — K2T/2 {Ki,K2 being 

constants), the line element can be written as 

ds^ = dT^ + h{z)dy^ + ^{dx + mdyf, (31) 

h[z) 2 

While the nonlinear electromagnetic field for < >> 1, from Eq. (p^), we see that 
v{t) goes to a constant. Then, from Eqs. ( [T^ ) for the fields D{t),B{t), we recover 
the Maxwell case. 



3.4. Relation to Einstein-Rosen Universes 

By means of a coordinate transformation, the metric (|l|) can be put in the form 
of an Einstein-Rosen line element Ii3 

ds^ = -XridZ^ - dT"^ + Gle-^dy'^ + (dx + mdyf]}. (32) 

The coordinate transformation for z Z depends on the value of e as follows 

z^{^~a)icoshZ-^, e = l, (33) 

z={!L + a)ismZ+^, e = -1, (34) 

e = 0, (35) 

We also note that for e = — 1, the functions in terms of the spatial coordinate z 
become periodic, then the topology of the universe in this case can be a closed one. 
The coordinate transformation for t, — > dT^, is not a simple one in the general 
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case due to the form of s{t), Eq. It involves elliptic integrals which lead to 

Jacobi family of elliptic functions cB and there is no analytical expression in terms 
of elementary functions. However, for particular cases such transformation can be 
given in a simple form. An example is for early times. In this case the relation 
between t and T is exp{^/eT) = t-'^(2et + cit^ + 2t^ e[e + cit + b'^t'^)). 

Another example is the case 6 = (vacuum) and e = 0. In this case the line 
element (js^) takes the form 

ds^ = ^{dZ"^ - dT^) + TZ[e-^dy^ + {dx + mdy f], (36) 

Gowdy constructed exact vacuum solutions of the Einstein field equations 
which represent inhomogeneous closed universes. These models possesses compact 
spacelike hypersurfaces as well as G2 invariance. The line element ( |36| ) can be writ- 
ten as a Gowdy model with thre-torus topology performing the coordinate trans- 
formation 

G = TZ^i, — + — (37) 
then one obtains the three-torus Gowdy line element 

ds^ = Zl— {de - de) + ^[e-'^dy^ + e^(dx + mdy)\ (38) 

If one arrives to a Gowdy model for vacuum, one can speculate about that a 
geometry so simple as Gowdy does not admit a Born-Infeld field. 

4. Concluding Remarks 

In this work we have presented a family of solutions to the Einstein-Born-Infeld 
equations for a space-time which is a G2 cosmological model. The spacetime de- 
scribes propagating gravitational plane waves coupled with a nonlinear spatially 
homogeneous e lectromagnetic field of the Born-Infeld type. We also present the 
limit in which the BI field become Maxwell electromagnetism. 

Some results are: 

The presence of an initial coordinate singularity depends on the choice of the 
conformal metric function Lp{t) = A cos{lt) + B sm{lt), for B 0, A = the solution 
exhibites initial singularity while ii B — Q, A there is no such coordi nate 
singularity. However, from the smooth behavior of the Weyl scalar, -02, we guess 
that there are no singularities at all. 

The presence of the BI field sets bounds on the amplitude of the periodic function 
(p{t). In this sense BI modifies the global expansion of the spacetime. If compared 
with the effect of linear electrodynamics, the presence of the BI field smooths the 
metric function s{t) (see Fig. 6). Since s{t) is involved in the expression of the 
scale expansion and of the Weyl scalars, then BI field smooths the curvature. The 
BI field also modifies quantitatively the curvature at i = 0, the effect being more 
sensible when the BI parameter b approaches the magnitude of a critic magnetic 
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field intensity H. We conjecture that a less restrictive spacetime (not type D) shall 
permit the existence of an inhomogeneous nonlinear electromagnetic field, whose 
consequences on the space curvature could be more drastic. 

For early times, t — > 0, the spacetime approaches a space of unpolarized grav- 
itational waves. For t >> 1 the spacetime becomes an inhomogeneous anisotropic 
spacetime with a homogeneous Maxwell field. 

For vacuum one obtains a Gowdy model of three-torus topology. 

It should be of interest to classify the spacetime in different regions according to 
the sign of the gradient of the function G(z, t) and see the distinct interpretations 
of the found solution according to each region in the {z,t) plane. 

Generalizing this solution to include a scalar field and to investigate the cou- 
pled effect with nonlinear electrodynamics for early universes could also lead to 
interesting results. 
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Fig. 2. It is displayed the metric function s{t) for different values of the BI constant b, b = 1, 1.5, 2, 3. For greater 
values of b the maximums are higher but the function preserves the shape of & > 2. These plots correspond to 
<p = 0.8sin(t), ci = 1 and e = 1. 
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Fig. 5. The function s{t) displays different behaviour at t = depending if Lp — 0.8sin(t) (continuous curve) or 
= 0.8cos(t) (dashed curve). 
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Fig. 6. The function s{t) for the hnear hmit corresponds to the dashed curve, the continuous hue is s{t) in the 
presence of BI field with 6=1. The constant ci = 1. s{t) of BI field is smoother than the one corresponding to 
linear electromagnetic field. The conformal metric function considered is 99 = 0.8sin(t). 



